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ABSTRACT 

This study reports an unusual heterogeneity in [^^C^®0]/[^^C^®0] abundance ratios of carbon 
monoxide observed in the gas phase toward seven ^ solar-mass YSOs and three dense foreground 
clouds in the nearby star-forming regions, Ophiuchus, Corona Australis, Orion, Vela and an isolated 
core, L43. Robust isotope ratios were derived using infrared absorption spectroscopy of the 4.7 /rm 
fundamental and 2.3 /xm overtone rovibrational bands of CO at very high resolution (A/AA « 95,000), 
observed with the CRIRES spectrograph on the Very Large Telescope. We find [^^C^®0]/[^^C^®0] 
values ranging from ^ 85 to 165, significantly higher than those of the local interstellar medium (~ 65 
to 69). These observations are evidence for isotopic heterogeneity in carbon reservoirs in solar-type 
YSO environments, and encourage the need for refined Galactic chemical evolution models to explain 
the ^^C/^^C discrepancy between the solar system and local ISM. The oxygen isotope ratios are con¬ 
sistent with isotopologue-specific photodissociation by CO self-shielding toward the disks, VV CrA N 
and HL Tau, further substantiating models predicting CO self-shielding on disk surfaces. However, 
we find that CO self-shielding is an unlikely general explanation for the high [^^C^®0]/[^^C^®0] ratios 
observed in this study. Comparison of the solid CO against gas-phase [^^C^®0]/[^^C^®0] suggests that 
interactions between CO ice and gas reservoirs need to be further investigated as at least a partial 
explanation for the unusually high [^^C^®0]/[^^C^®0] observed. 

Subject headings: astrochemistry — protoplanetary disks — stars: solar-type — stars: protostars — 
infrared: stars — molecular processes 


1. INTRODUCTION 

Observations of carbon monoxide (CO) provide a ro¬ 
bust window into circumstellar chemistry. Due to the 
large percentage of carbon (~ 15 to 30%) and oxygen 
10%) tied up in CO in the early stages o f young 
stellar object (YSO) formation (Dickman |1978 ), CO is 
a valuable tracer of chemical inheritance from the enve¬ 
lope to the disk in solar-type YSOs. The use of high- 
resolution spectroscopy to observe CO at infrared wave¬ 
lengths enables high-precision investigations of isotope 
fractionation in both oxygen and carbon toward YSOs 
in embedded and disk-phase evolutionary stages, and is 
consequently an important tool for connecting the so¬ 
lar system to currently forming protoplanetary systems. 
Examples include observational studies toward diffuse 
clouds and inferring isotopologue-specific photodissoci¬ 
ation by CO self-shielding as a potentially important 
mechanism of carbo n isotope fractionation in CO in the 
inters tellar medium (Bally & Langer||1982 Sheffer et al. 
2002|). Further, high-resolution near-infrared observa- 


tions of th e CO rovibrational band toward protop lane- 
tary disks (Brittain et al. 2005 Smith et al. 2009) and 
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protostellar envelopes have recently revealed signatures 
in CO isotopologue ratios that are consistent with CO 
self-shielding, with precisi on sufhciently high for com¬ 
parison to meteoritic data (Smith et al. 2009). 


Observations of isotope ratios toward Y SOs and, where 
possible, their parent clouds, are valuable in assessing 
theorized similarities between circumstellar envelopes, 
disks and dense clouds local to their origin. This is im¬ 
portant not only for evaluating predictions of galactic 
chemical evolution (GCE) models that are constrained 
by solar system abundances, but also for characteriz¬ 
ing the range of carbon isotope ratios toward YSOs 
varying with evolutionary stage, and establishing ac¬ 
curate astrochemical constraints on protoplanetary en¬ 
vironments. Observational studies can help demon¬ 
strate whether there are significant isotopic differences 
or trends between diffuse clouds, dense molecular clouds, 
and YSOs. These in turn help construct explanations for 
apparent chemical anomalies, such as the ^^C/^^C solar 
system-interstellar medium (ISM) discrepancy and solar 
system oxygen isotope anomaly. 

This paper focuses on the carbon isotopic systems 
in YSOs as traced by high-resolution observations of 
CO isotopologues, with additional mention of significant 
findings in oxygen isotopes. Current measurements of 
the ^^C/^^C ratio for the solar system are ~ 86 to 89, de¬ 
rived from photospheric mo dels applied to the solar CO 


rials ( 

Rosman &: 'i'aylor 

1998 

Wilson|1999l 

Coplen et al. 

to 

o 

o 

to 

Clayton & Wittier 

2U04 

Young et al.| 

Oil). Values 


of ^^C/^'®C for the solar system have been consistently 
higher than those of ^ 62 to 69 observed toward the lo¬ 
cal ISM, measured using millimeter to sub-millimeter ro¬ 
tational transitions of various carbon-bearing molecules, 
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CO, CN, CH+, and H2CO ( L anger fc Penzias|l993 Wil- 
son||1999 Milam et al.||2005 ' Casassus et al.||20U5| )' 'I’he 


apparent discrepancy in between the solar sys¬ 

tem and ISM has been attributed to GCE leading to 
enrichment of relative to over the ~ 4.6 Gyr life- 


TABLE 1 

Summary of Observed Targets 


time of the solar system (e.g. Milam et al.|2005 Prantzos 


et al. 1996). An important property of GCE models is 
the assumption that the solar system records the chemi¬ 
cal state of the molecular cloud material ^ 4.6 Gyr ago. 
That is, the assumption that the solar nebula was iso- 
topically representative of its parent cloud at its Galac¬ 
tic radius at the time of formation should also hold for 
present-day solar-type YSOs. 

The challenge for observational comparisons of so¬ 
lar system isotope ratios in general is the need for 
high precision in molecular abundance measurements. 
While techniques using strong rotational lines in the 
millimeter/sub-millimeter range often involve model- 
dependent comparisons between optically thick and thin 
lines, high-resolution infrared spectrographs have en¬ 
abled significant breakthroughs in observations of nearby 
YSOs. Infrared absorption spectroscopy of optically thin 
CO lines allows for precise column density measurements 
along pencil-beam lines-of-sight, enabling the derivation 
of precise column densities and isotope ratios from the 
line profiles of as many as four currently observable CO 
isotopologues (^^Ci^O, i2c18o, 12^170) 

bust isotope ratios so derived can be compared with sta¬ 
tistical significance to sub-millimeter ISM values, solar 
abundances, and meteoritic oxygen. A high degree of 
precision (within 10%) on the absolute values of the col¬ 
umn densities has been dem onstrated for the fou r most 
common CO isotopologues (Smith et al. 2009). The 
study reported here builds on this work. 

Our targets span a range in ^20160 optical depth, 
from nearly saturated to undetected, providing an oppor¬ 
tunity to investigate a potential observational connection 
between the CO ice and gas reservoirs by comparing the 
gas-phase [^^C^®0]/[^^C^®0] to the CO ice column along 
a single line-of-sight. Since solids tend to concentrate the 
heavy isotope, relative ice-gas partioning may lead to a 
paucity of ^^C (or excess of ^^C) in the gas phase, with 
the opposite signature sequestered in the ice. 

2. OBSERVATIONS 

Data for all targets were obtained as part of a Euro¬ 
pean Southern O bservatory Large Program (Pontoppi- 
dan et alT]|2011b ) to observe ~ 100 YSOs and protoplan¬ 
etary disks with the Cryogenic Infrared Echelle Spectro¬ 
graph (CRIRES) at the Very Large Telescope (VLT) in 
Chile. The CRIRES instrument is an adaptive optics 
(AO)-assisted spectrometer that operates at very high 
resolving power (A/AA « 95,000). 

2.1. Strategy 


Target 

Region^ 

Distance 

(pc) 

Comments^ 

DoAr24E S,N 

Oph 

125 

Stage II/IR companion 

IRS 43 

Oph 

125 

Stage I 

IRS 51 

Oph 

125 

Extincted stage II 

IRS 638 

Oph 

125 

Stage I 

RNO 91 

L43 

125 

Stage I 

VV CrA N,S 

CrA 

130 

Stage II/IR companion 

HL Tau 

Tau 

140 

Stage I 

Reipurth 50 

Ori 

470 

Stage I/FU Ori 

LLN 19 

Vel 

700 

Stage I 


Note. — 

The evolutionary st age refers to the system proposed by |Ro-| 
|bitaille et al. 1 1|2006||. The two targets with N/S components 
are binaries, tor w hich the secondary i s a highly extincted 
“infrared companion” ([Koresko et al.|1997[|. 


8 IRS 63 may be formed from an isolated core jRidge et al. 

[20^ . - 

of the fundamental band by a factor of ~ 100, the for¬ 
mer are typically optically thin for thereby en¬ 

abling the precise measurement of column densities for 
this molecule. The spectra of the fundamental rovibra- 
tional bands were obtained using the 0"2 slit, resulting in 
a resolving power of i? = A/AA « 95, 000 (corresponding 
to 3.2kms“^). Several settings in the fundamental band 
were obtained for all targets in order to probe a range of 
rotational levels for ^20180 I20i70^ span¬ 

ning J = 0 to at least J = 9. Observations were timed 
such that the Doppler shifts of the targets relative to 
the Earth’s atmosphere were at least 20kms“^ due to 
a combination of the intrinsic target velocity and reflex 
motion of the Earth. As a consequence, telluric CO lines 
are typically not affecting the absorption lines of the tar- 
gets. 

Table summarizes the region, distance and evolu¬ 
tionary stage of the observed targets, listed in order of 
increasing distance from the local solar neighborhood 
(i?GC 8 kpc). Specific observing parameters are listed 
in Table Selected targets probe a range of circumstel- 
lar environments, as noted. Because the lines-of-sight 
traced by infrared absorption may pass through the up¬ 
per layers of an inclined disk as well as through envelope 
material and a foreground molecular cloud unrelated to 
the YSO, some portion of each of these environments may 
be represented in the integrated spectra, depending on 
the target. For instance, near-infrared absorption spec¬ 
tra observed toward the extincted protoplanetary disk, 
IRS 51, will be an integrated sight line consisting mostly 
of the foreground Ophiuchus ridge and a relatively small 
contribution of the protoplanetary disk. This geom etry 
has been predicted by interferometric observations (van 


Kempen et al.|2009 ) and is supported by our observations 
that reveal cold-only gas lines. 


Targets tracing high column densities of molecular gas 
were selected based on evidence of deep GO absorption 


lines at 4.7 ym in si 

lectra obtained with VLT-ISAAC 

(IPontoppidan et al.| 

2003) and Keck-NIRSPEG (e.g., 

Boogert et al. 

20U2b 

). belected targets were observed 


by CRIRES in the GO overtone (u = 2 — 0) band at 2.3 
fim and the fundamental (u = 1 — 0) band at 4.7 
Because the overtone transitions are weaker than those 


2.2. Data reduction 

The CRIRES spectra were reduced using standard 
procedures for infrared echelle spectroscopy, including 
flat-held correction, adjustments to account for detec¬ 
tor non-linearity, and linearization of the spectral traces 
in dispersion- and cross-dispersion directions. The spec¬ 
tra were optimally extracted, wavelength-calibrated us¬ 
ing the telluric absorption lines referenced to an atmo- 
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TABLE 2 

Summary of Observed Targets 


Target 

RA 

{J2000) 

Dec 

(J2000) 

Spectral Range 
(Xtm) 

Obs. Date 
(UT) 

DoAr24E S,N‘‘ 

16 26 23.36 

-24 20 59.8 

4.644 - 
2.338- 

4.768 

2.396 

2007-09-02 

2010-03-18 

IRS 43 

16 27 26.94 

-24 40 50.8 

4.645 - 
2.338 - 

4.900 

2.396 

2008-08-05 

2008-08-04 

IRS 51 

16 27 39.83 

-24 43 15.1 

4.645 - 
2.338 - 

4.768 

2.391 

2008-08-02 

2008-08-03 

IRS 63“* 

16 31 35.65 

-24 01 29.3 

4.645 - 
2.338- 

4.978 

2.396 

2007- 04-25 

2008- 08-07 

RNO 91 

16 34 29.28 

-15 47 01.5 

4.645 - 
2.338- 

4.855 

2.396 

2010-03-18 

2010-04-08 

VV CrA N,S 

19 03 06.74 

-37 12 49.7 

4.701 - 
2.338 - 

4.815 

2.396 

2007-08-31 

2007-08-31 

HL Tau 

04 31 38.44 

18 13 57.7 

4.643 - 
2.337- 

4.853 

2.396 

2010-01-30 

2010-01-28 

Reipurth 50 

05 40 27.70 

-07 27 28.0 

4.645 - 
2.338- 

4.901 

2.396 

2007-10-11 

2007-10-17 

LLN 19 

08 48 48.00 

-43 32 25.0 

4.644 - 
2.253 - 

4.867 

2.355 

2007-04-25 

2007-04-25 


Note. — 

^ The non-coevality of these observations may contribute some level of un¬ 
certainty to the derived isotopic abundances, as discussed in the text. 




Fig. 1. — Infrared absorption spectra of the CO fundamental 
rovibrational band toward the DoAr24E N,S binary and the disk 
VV CrA N, which have weak CO ice optical depth upper limits 
(Tabled. Spectra for DoAr24E N,S show few CO absorption lines, 
while VV CrA N shows a rich CO forest of narrow rovibrational 
absorption lines. Sample lines are marked. The low-level emission 
from hot CO gas toward DoAr24 S and VV CrA N is possibly 
representative of inner disk disk material. The weak emission line 
at ~ 4.65 fim for DoAR24E N is HI Pfund/3. The main portions of 
the spectra used for the line analyses are presented. 


spheric model spectrum, and transformed to the local 
standard-of-rest frame. Relative flux calibration was car¬ 
ried out by dividing the targets by observations of early- 
type telluric standards observed immediately surround¬ 
ing observations of the science target at the same air- 
mass. The division was optimized by modifying the tel¬ 
luric spectra to account for any difference in grating an¬ 
gle, airmass, and spectral resolving power. To aid in the 
measurement of any intrinsic line broadening, the resolv¬ 
ing powers for each science spectrum were measured from 
telluric spectra to within an accuracy of ~ 0.1 kms“^, by 
fitting the “Reference Forward Model” (RFAiQ see also 
Section 3.2). F urther details on CRIRES data pro cessing 
is described in Pontoppidan et al. (2008[ [^lla[). 


3. METHODS AND RESULTS 
3.1. CO spectra 

The fundamental (u = 1 — 0) and first-overtone {v = 
2 — 0) rovibrational bands are shown for all targets in 
Figures |T] through in ascending order of CO ice opti¬ 
cal depths at ^ 4.675/xm. A forest of narrow CO rovi¬ 
brational absorption lines is seen to varying degrees in 
all spectra. A magnified spectral region of HL Tau is 
shown as a rraresentative example of individual line pro¬ 
files (Figurej^. Broad, complex {v = 1 — 0) CO emission 
is often seen in the spectra, and likely ori ginates in the 
inner most (i? < 1AU) regions of the disks |Herczeg et al.j 
20111. However, the emission does not interfere with the 
analysis, as it is much broader than the absorption, and 
is generally weak or absent in the rare isotopologue tran¬ 
sitions. The potentially affected (v = 1 —0) ^^C^®0 lines 
are not used for column density measurements. Broad 
absorption observed toward several objects at (u = 2 —0) 
(Figures [^1^ likely originates in the stellar photosphere 
and is clearly separated from the narrow interstellar CO 
lines. 

Temporal variation in [^^CO]/[^^CO] was found toward 
VV CrA S, with a [^^CO]/[^^CO] abundance ratio of 


® http://www.atm.ox.ac.uk/RFM 
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Fig. 2. — Infrared absorption spectra of the CO fundamental 
rovibrational band toward VV CrA S, HL Tau and LLN 19, char¬ 
acterized by weak lower-limits fVV CrA S and HL Tau), or shallow 
(LLN 19) CO ice depths (TableQ. HL Tau and LLN 19 have a rich 
CO forest of narrow rovibrational absorption lines due to the vari¬ 
ous CO isotopologues present. Sample lines are marked. Few lines 
of the rarer isotopologues (^^C^^O and ^^C^^O) are seen toward 
the more deeply embedded object, LLN 19. Low-level emission to¬ 
ward VV CrA S and HL Tau are possibly representative of hot CO 
gas of the inner disk. 


136 ± 1.3 when using a (v = 1 — 0) observation four 
months prior to the first-overtone (u = 2 — 0), as com¬ 
pared to 127 ± 1 when coeval observations were used (er¬ 
rors are statistical la values). The possibility for tempo¬ 
ral spectral variability should therefore be considered for 
DoAr24E and IRS 63, targets for which only non-coeval 
observations are available. 


3.2. Derivation of column densities and isotopologue 

ratios 

Total abundances for each isotopologue were derived 
by measuring the column density of CO in each observed 
rovibrational level of the ground vibrational state which, 
in turn, was obtained from the intrinsic optical depths 
at line center (to) derived from a best-fit absorption line 
profile to each rovibrational transition. Individual col¬ 
umn densities were used to make a best fit to a Boltz¬ 
mann population, then summe d over all levels. This 


method is described in detail in Smith et al. (2009). 

The technique of using the line prohles to directly de¬ 
termine column densities takes advantage of the high 
spectral resolution of CRIRES to spectrally resolve the 
velocity dispersion of the individual CO lines. This in 



Fig. 3. — Infrared absorption spectra of the CO fundamental 
rovibrational band toward RNO 91 and Reipurth 50, with moder¬ 
ate optical depths of CO ice (r <l).The rich CO forest of narrow 
rovibrational absorption lines are due to the various CO isotopo¬ 
logues present. Sample lines are marked. RNO 91 shows low-level 
emission possibly representative of the inner disk. 



Fig. 4. — Infrared absorption spectra of the CO fundamental 
rovibrational band toward IRS 63, IRS 43, and IRS 51, embedded 
YSOs with deep CO ice features. The CO absorption toward IRS 
51 is particularly cold, leading to detection of only low-J transi¬ 
tions. Sample lines are marked. 

turn leads to greater precision than using unresolved lines 
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2.34 2.35 2.36 2.37 2.3B 2.39 2.335 2.340 2.345 2.350 2.355 

Wavelength (pm) Wavelength (pm) 


Fig. 5.— Infrared absorption spectra of the CO rovibrational first-overtone bands toward the YSOs shown in Figures ^ andSeveral 
detected rovibrational lines are marked. Broad absorption lines at ^ 2.34 — 2.353 pm, most evident toward DoWR24E N,S, VV 

CrA S, and LLN 19, are due to photospheric CO. 


which require a curve-of-growth analysis to derive intrin- 
sic line widths dMitchell et al. 1989[ 19^ 1991| Sheffer et 


at the time of observation. The atmospheric molecular 

__ _ _ parameters were taken fr om the HITRAN 2008 database 

^|1992 Mitchell & Maillard 1998 ' Lanibert et al.||1994[ (Rothman et al. 2009). There are small relative dif- 


Shefter et al.| 120021 [Federman et ^ |2003[ [Goto et al.| 

2003| |Brittain et al'||2005p. With CRIRES, the resolving 
power is suthciently high to partially resolve the intrinsic 
line profiles (FWHMintrinsic ~ FWHMinstrumentai): albeit 
marginally in some cases: VV CrA S, RNO 91, and IRS 
43. Therefore, it is possible that the derived column den¬ 
sities toward these targets may be more uncertain. 

The measured instrumental broadening, intrinsic ve¬ 
locity dispersions and optical depth ranges for each ob¬ 
served isotopologue are shown in Table Line-by-line 
optical depths and Doppler shifts for observed targets are 
available in the online version of this paper, with a sam¬ 
ple partial line list shown in Table Complete data ta¬ 
bles are available in the online version of the papeiFj In¬ 
strumental broadening values were modeled individually 
for each spectrum in its entirety, following the method in 
Seifahrt et al. (20101 for using an atmospheric model for 
telluric spectra. The Reference Forward Model (RFM; 
http://www.atm.ox.ac.uk/RFM/) code was used to de¬ 
rive a detailed fit to the individual telluric standard star 
spectra, with consideration of meteorological conditions 

^ Complete data tables can be temporarily found at 
http://naturalsciences.org/data/apj2015 


ferences in resolvihg power between the standard stars 
and science targets due to differences in seeing or adap¬ 
tive optics (AO) correction, but these differences were 
accurately measured and corrected for by minimizing the 
continuum residuals after telluric division. 

Optimized intrinsic velocity dispersions were derived 
from the weighted mean of dispersions measured in indi¬ 
vidual ^^C^®0 lines, assuming that both the intrinsic and 
instrumental profiles are well-approximated by Gaussian 
profiles. The ^^C^®0 lines were used because this rare 
isotopologue is always optically thin, and therefore is not 
subjected to saturation line broadening, which may affect 
the more common ^^CO isotopologue. Optical depths of 
the optically-thin (u = 2 — 0) lines were sim¬ 

ilarly computed using the derived intrinsic broadening 
from the (w = 1 — 0) lines. It can be noted 

that, even when the resolving power of the two spectral 
bands differed for a single target, we obtained good fits 
to lines in both bands under the reasonable assumption 
that the intrinsic profile width is independent of transi¬ 
tion and isotopologue. Further, it is reasonable to use 
the same velocity dispersion parameter for a target with 
a two-temperature distribution of the data, since thermal 
contributions are very small contributors to the total in- 
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2,340 2.345 2,350 2.355 2.360 

Wavelength (\im) 


Fig. 6. — Infrared absorption spectra of the CO rovibrational first-overtone bands toward the YSOs shown in Figuresand [4| Several 
detected rovibrational lines are marked. Broad photospheric CO lines are observed toward RNO 91 and more wSkly toward IRS 

63. 



Wavelength (pm) 

Fig. 7. — Representative magnification for fundamental (top) and 
first-overtone (bottom) bands, here shown for HL Tau. Isotopo- 
logue lines are marked. 

trinsic dispersion (0.077 and 0.242 km s“^ for at 

20 and 200 K, respectively). 

The rotational analyses for all targets are shown in Fig¬ 


ures and Sub-level column densities, Nj/{2J + 1), 
derived from the fits to the absorption lines, are plot¬ 
ted against the energy, E/k^ of the rotational transition 
state, J. For targets where cold gas dominates, a single¬ 
temperature model is sufficient to fit the data, while tar¬ 
gets with absorption that includes higher-J transitions 
are best fit with a two-temperature model, representing 
cold and warm environments along the line-of-sight, and 
allowing for different excitation temperatures for each 
isotopologue along a given line-of-sight. Derived isotopo- 
logue temperatures from the fits are shown in Table 
High temperatures (~ 80 to 250 K) are roughly consis¬ 
tent with those of the inner parts of typical YSO en¬ 
velopes, while low temperatures (~ 5 to 30 K) may trace 
outer envelopes or molecular cloud material in front of 
the YSOs. For lines-of-sight with no warm component, 
it is likely that the gas is from a foreground cloud not 
related to the YSO itself (DoAr24E N and S), VV CrA 
S, and IRS 51. 

Total column densities and isotopologue ratios are 
shown in Table Due to the dominance of higher-J lines 
in targets witliboth low- and high-temperature lines, 
the higher-temperature ratios are generally be tter con¬ 


strained within two-temperature environments (Smith et 


al. 


2009). 


Statistical errors on the column densities are gener¬ 
ally found to be small, reflecting the very high signal- 
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TABLE 3 

Summary of observed line parameters 


Target 

FWHMinstrumental^ 

M (km s“^) K (kms“^) 

FWHMintrinsic® 

(kms 

t ( 12 c 160 ) 

v = (2-0) 

t ( 13 c 1 ® 0 ) 

•u = (1 - 0 ) 

t ( 12 c 180 ) 

i; = (1 - 0 ) 

r; = (1 - 0 ) 

DoAr24E S 

2.79 

3.12 

2.33 ±0.56 

0.184 - 0.353 

0.217- 1.033 

0.027- 0.077 

- 

DoAr24E N 

2.79 

3.12 

2.33 ±0.56'= 

0.075 - 0.330 

0.198- 1.014 

0.047 - 0.068 

- 

VV CrA N 

2.82 

3.49 

4.28 ±0.07 

0.064 - 0.573 

0.568- 1.322 

0.095 - 0.174 

0.020 - 0.038 

VV CrA S 

2.82 

3.49 

1.74 ±0.12 

0.170 - 0.614 

0.225 - 1.377 

0.119 - 0.228 

0.021 - 0.034 

HL Tau 

3.26 

3.43 

2.63 ±0.05 

0.027-0.987 

0.036 - 2.896 

0.079 - 0.310 

0.021 - 0.052 

LLN 19 

3.22 

4.68 

3.52 ±0.27 

0.196 - 0.437 

0.627-0.909 

0.072 - 0.108 

- 

RNO 91 

3.26 

3.72 

1.67 ± 0.27 

0.105 - 0.843 

0.089- 2.290 

0.046 - 0.370 

0.059 - 0.072 

Reipurth 50 

3.44 

4.54 

3.08 ±0.02 

0.225 - 1.366 

0.338- 2.831 

0.027-0.707 

0.044-0.138 

IRS 63 

3.57 

3.39 

3.83 ±0.35 

0.057-0.811 

0.171 - 1.582 

0.141 - 0.245 

0.043 - 0.081 

IRS 43 

3.03 

3.45 

1.86 ±0.11 

1.076 - 2.563 

0.212 - 4.542 

0.073 - 1.115 

0.033 - 0.280 

IRS 51 

3.03 

3.37 

1 .88'^ 

0.504 - 0.822 

0.739- 1.437 

- 

- 


Note. — 

^ The instrumental broadening is measured by fitting an atmospheric model to unresolved telluric lines. The uncertainty is 
< 0.1 km s“^. These refined mea surements of instru mental broadening account for the slight variation in results for VV CrA 
and Reipurth 50 as compared to ( [Smith et al.|200^ . 

^ The intrinsic velocity dispersion is the Gaussian FWHM before being convolved with the measured instrumental CRIRES 
broadening for the individual spectra. It is based on a weighted mean FWHM from a selection of lines, unless otherwise 

noted. Uncertainties reflect the weighted means of the individual line widths. 

^ Mean value from DoAr24E S lines R(2) and R{1). 

^ This width is estimated based and lines, since the source has very weak lines. 


TABLE 4 

Optical depths (to) and Doppler shifts for CO 

ISOTOPOLOGUE ROVIBRATIONAL LINES USED IN DERIVING 
COLUMN DENSITIES FOR DoAr24E 


Isotopologue Transition 

ro 

y^sR 

(km s~^) 

12C160 

(2,0) R(3) 

0.279 ± 0.002 

3.09 ±0.01 

12 c 16 o 

(2,0) R{2) 

0.353 ±0.002 

3.31 ± 0.01 

12 c 16 o 

(2,0) R{1) 

0.326 ±0.002 

3.42 ±0.01 

12 c 16 o 

(2,0) R(0) 

0.262 ±0.002 

3.84 ± 0.02 

12 c 16 o 

(2,0) P(l) 

0.184 ± 0.002 

3.53 ±0.02 

12C160 

(2,0) P(2) 

0.255 ±0.002 

3.47 ± 0.01 

12C160 

(2,0) P(3) 

0.250 ± 0.003 

4.06 ±0.03 


Note. — 

^ This is a partial table for DoAr24E S. Line-by-line ta¬ 
bles for all observed targets are available in the online 
version of the paper, and can be temporarily found at 
http://naturalsciences.org/data/apj 2015 
^ Uncertainties are Icr derived from the line fits. 


to-noise ratio of the data. For example, average signal- 
to-noise per pixel in the fundamental bands is ^ 250. 
With an average CRIRES FWHM of 2.0 pixels, we re¬ 
trieve average line errors of ^ 0.3% (1/(250 x -s/^)). 
First-overtone bands average ~ 1%. We therefore expect 
to know the column densities to a fraction of a percent, 
which is represented in our results. These errors do not 
take into account any potential systematic errors in ex¬ 
tracting the line strengths and profiles due to inherent 
non-gaussianity of the line-fitting procedure, which we 
cannot quantify with certainty. However, to roughly es¬ 
timate the magnitude of systematic retrieval errors, we 
carried out a small number of forward models with 2- 
dimensional radiative transfer models (RADLite 


Pon- 


toppidan et al. 2009) simulating synthetic CO isotopo- 
logue absorption spect ra us i ng a basic one-dimensional 
infalling envelope from Shu] (1977). The synthetic spec¬ 
tra were convolved with the appropriate instrumental 
resolution and isotopologue ratios retrieved from the syn¬ 
thetic spectra using the same methods as for the observed 


data. In each case, we were able to reproduce the input 
ratios to better than 10%. We therefore estimate that 
model-dependent systematics add at most a 10% error 
to the retrieved isotopologue ratios, in addition to the 
statistical error. 


3.3. Notes on individual sources 

The DoAr24E system is a hierarchical multiple for 
which we have spectra of two components, with the ro¬ 
tational analyses nearly identical for both components. 
This lack of signihcant spatial variability strongly sug¬ 
gests that we are tracing spatially extensive and homoge¬ 
nous foreground material, that in turn is the dominant 
contributor smoothly distributed on scales of at least a 
few 100 AU. The lack of isotopic variability between the 
DoAr24E components gives us further confidence on our 
retrieval method for deriving precise isotope ratios. 

Due to high extinction leading to a faint v = {2 — 
0) band continuum for the embedded YSO IRS 43, we 
detect high-J lines for all CO isotopologues in = (1—0), 
but only low-J lines for ^^C^®0 in the r; = (2 — 0) band. 
Given this disparate set of temperature ranges, our data 
for IRS 43 are therefore not useful for measuring CO 
isotopologue abundance ratios relative to 

Toward VV CrA S, we derive a ^^C^®0 column density 
of 1.178±0.004x 10^® cm“^. This is nearly three orders of 
magnitud e higher than the val ue of 2.6± 0.4 x 10^® cm“^ 
derived in Kruger et al. (2011). Our rotational tempera¬ 
ture (14 ± 1K) i s also signihcantly lo wer than their value 
(29±4K). Since |Kruger et al.|(|2dll|) used the (r; = 1 —0) 
spectra from the archive, these dinerences could be ex¬ 
plained by strong saturation effects resulting from their 
use of the highly optically thick u = (1 — 0) ^^C^®0 
lines. This difference reinforces the necessity of using 
V = (2 — 0) data for deriving precise ^^C^®0 column den¬ 
sities. 


Smith et al. (2009) proposed that the cold gas compo- 
nent observed toward VV CrA S may be interpreted as 
either molecular cloud material in front of the VV CrA 
binary system, or disk gas toward the primary compo- 
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TABLE 5 

Measured CO isotopologue rotational temperatures 


Target 

Component^ 

Texfi^cieo) 

(K) 

Tex(13cl60) 

(K) 

7^ 

X 

00 

o 

Texfl^clTO) 

(K) 

DoAr24E S 

single 

30 ± 1 

16 ± 1 

8 ± 2 

- 

DoAr24E N 

single 

29 ± 1 

17 ± 1 

13 ± 12 

- 

VV CrA N 

high 

249 ± 2 

255 ±8 

179 ±5 

207 ± 39 

VV CrA N 

low 

14 ± 1 

14 ± 1 

6 ± 1 

13 ±3 

VV CrA S 

single 

14 ± 1 

12 ± 1 

8 ± 1 

9± 1 

HL Tau 

high 

91 ± 1 

97 ± 1 

89 ±6 

109 ± 17 

HL Tau 

low 

24 ± 1 

18 ± 1 

25 ±3 

19 ±5 

LLN 19 

single 

28 ± 1 

15 ± 1 

23 ±2 

- 

RNO 91 

high 

95 ± 1 

100 ± 1 

121 ± 21 

- 

RNO 91 

low 

19 ± 1 

13 ± 1 

29 ±5 

37 ± 6 

Reipurth 50 

high 

148 ± 8 

204 ± 1 

149 ± 1 

112 ±3 

Reipurth 50 

low 

20 ±3 

16 ± 1 

21 ± 1 

14 ± 1 

IRS 63 

high 

145 ±3 

144 ± 4 

91 ± 8 

76 ±40 

IRS 63 

low 

22 ± 1 

17 ± 1 

- 

- 

IRS 43 

high 

27 ± 1 

82 ± 2 

106 ± 15 

88 ±21 

IRS 43 

low 

- 

8 ± 1 

13 ± 1 

12 ± 2 

IRS 51 

single 

23 ± 1 

16 ± 1 

- 

- 


Note. — 

All uncertainties are Icr statistical error derived from the rotational analysis. 

^ The temperature component. Single if a single temperature component was used. High 
and low for two-component fits. 


nent as suggested by Kruger et al. (20111. We point out 
that the former is niore likely because the presence of 
both warm and cold components are expected for a disk, 
while a foreground cloud will typically only show a single 
cold component, such as that seen towards VV CrA N. 
Taken together, the likely geometry for the VV CrA sys¬ 
tem and the spectroscopic results render it unlikely that 
the VV CrA S spectra represent outer disk. VV CrA 
N observations most likely represent warm and cold disk 
gas. 

3 . 4 . Carbon isotope ratios versus Galactocentric radius 
Derived [^^C^®0]/[^^C^®0] ratios a re p lotted against 


Galactocentric radius (Rgc) in Fi gure 10 The solar sys¬ 
tem a nd local ISM values of ~ 87 (Scott et al.|2006|) and 
^ 68 (Milam et al. 2005) are plotted at 8 kpc, respec¬ 
tively. Derived [^^C^'^0]/[^^C^®0] toward our targets 
range from ~ 85 to 165, values inconsistent with those 
toward most YS Os and molecular clouds, or from GCE 
dPrantzos et al. 


1996) and regression models (Milam et 


2(105 r 'I’he distribution of data from the literature 
approximates the general trend of increasing metallicity 
toward the Galactic center, with the corresponding in¬ 
crease in ^^C/^^C with greater Rgc- 

3 . 5 . Derivation of CO ice fractions 

Previous studies have discussed the exchange between 
CO solid- and gas-phase reservoirs in disks and en velopes 
(e.g.|van Dishoeck et al.|1996||Boogert et al |2000|). With 
a range ot CO ice column densities centered at 4.75 ^m 
(Figures 0 through 1^, our data set provided an oppor¬ 
tunity to investigate relationships between the gas-phase 
isotopologue ratios and the CO ice fraction along single 
lines-of-sight. CO ice column densities were calculated 
fro m optical depth an alyses in Pontoppidan et al. (2003) 
and 
O 


Thi et al. (2010). 

: the broad CO ice profile, the pure, cold CO compo¬ 


nent is most relevant with respect to interactions with 
CO gas, as it is the solid state reservoir that would 
most likely directly exchange with the gas in cold, dense 


clouds. We used optical depths for this component to de¬ 
rive total CO ice column densities for each targe t. Using 
the phenomenologi cal decomposition analyses in Pontop-| 
pidan et al. (2003), the pure component correspon ds to 
the middle optical depth (rmc) profile measurement; [Thi| 
et al.| (2010) was used for LLN 19. For targets with no 


(VV CrA, DoAr24E) or very weak (HL Tau) spectral ice 
features, we used upper limits from the CRIRES funda¬ 
mental bands. 

Total column densities of pure CO in the ice phase 
(Vpiirs noF,) we r e com puted using the relation from Pon¬ 
toppidan et al. (2003), 


-^Pure.CDE = 6.03 Cm ^ X T, 


^ ^bulk> 


( 1 ) 


where Tmax is the optical depth at line center of the pure 
CO component, Abuik is the band strength of ^^CO, and 
the numerical factor, 6.03 cm“^, considers the effect of 
grain shape on the band strength. In this study, we used 
the continuous distribut ion of ellipsoids (CDE) model 
( Pontoppidan et al.|2003 ) in computing solid CO column 
densities. 

Optical depths, column densities, and total CO ice 
fractions are listed in Table in ascending order of CO 
ice optical depth. Although it is assumed that cold CO 
gas will more likely exchange with CO ice, for targets 
where CO gas was observed in both cold and warm en¬ 
vironments, values for CO ice fractions are shown sepa¬ 
rately for each environment. 

The relation between the CO ice fraction and 
j12c16o]/[13c16o] gas-phase ratios is shown in Fig¬ 
ure with CO ice fractions plotted against derived 
[i2CfD]7[i3CO] for cold (T ~ 5 to 30 K) and warm 
{T ~ 80 to 250 K) regimes, the latter only for targets 
where warm gas was observed. We distinguish between 
targets where only cold CO gas is observed, and targets 
where both warm and cold temperatures is found. In the 
latter cases, cold-gas measurements are likely more un¬ 
certain due to the contribution of the high-temperature 
on the low-temperature lines. 
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TABLE 6 

Derived Isotopologue Column Densities and Isotopologue Ratios 


Target 

Comp. 

jV(12c16o) 
(IQl® cm-3) 

[12c16o]/[13c160] 

[I2ci®0]/[i2ci®0] 

[i2ci®o]/[i2ci'^o] 

[12c18o]/[12c170] 

DoAr24E S 

single 

1.098 ±0.005 

98 ± 1 

:2500 ± 560 

- 

- 

DoAr24E N 

single 

1.094 ± 0.003 

102 ± 3 

:2200 ± 1620 

- 

- 

VV CrA N 

high 

7.618 ±0.074 

95 ± 1 

740 ± 15 

3085 ± 440 

4.1 ± 0.6 

VV CrA N 

low 

1.413 ±0.023 

127 ± 2 

1750 ± 70 

4000 ± 630 

2.3 ±0.4 

VV CrA S 

single 

1.178 ±0.004 

127 ± 1 

:1150±20 

:3850 ± 580 

3.4 ±0.5 

HL Tau 

high 

4.467 ± 0.064 

85 ± 1 

795 ± 80 

3500 ± 660 

4.4 ±0.9 

HL Tau 

low 

1.067 ± 0.016 

74 ± 1 

760 ± 80 

2670 ± 370 

3.5 ± 0.6 

LLN 19 

single 

2.020 ± 0.034 

137 ± 3 

:1150±90 

- 

- 

RNO 91=" 

high 

2.364 ± 0.032 

92 ± 2 

740 ± 210 

2170 ± 280 

2.9 ±0.9 

RNO 91=" 

low 

0.689 ±0.007 

80 ±2 

350 ± 40 

630 ± 80 

1.8 ±0.3 

Reipurth 50 

high 

10.81 ± 0.766 

124 ± 9 

580 ± 40 

2480 ± 200 

4.3 ±0.2 

Reipurth 50 

low 

1.976 ±0.157 

76 ±6 

420 ± 30 

3470 ± 380 

8.2 ± 0.6 

IRS 63 

high 

5.743 ±0.151 

145 ± 6 

690 ± 60 

1830 ±510 

2.6 ±0.8 

IRS 63 

low 

1.904 ± 0.026 

110 ± 3 

- 

- 

- 

IRS 43'" 

high 

6.908 ±0.158 

:210± 10 

:1160± 180 

:2600 ± 960 

2.2 ± 0.9 

IRS 43'" 

low 

- 

:300± 10 

:1200±40 

:4500 ± 580 

:3.7± 0.5 

IRS 51 

single 

2.026 ±0.036 

165 ± 10 

- 

- 

- 

Local ISM= 

- 

- 

69 ±6 

557 ± 30 

2005 ± 155 

3.6 ±0.2 


Note. — 

All uncertainties are Icr statistical errors derived from the line fits. 

Isotope ratios that may be affected by significant systematic uncertainty are indicated with a . 

^ The lines for RNO 91 were fit with a single-temperature model. 

Ratios relative to are uncertain because only low-J lines are available, while the low-J lines for the other 

isotopologues are poorly con strained. 

Values from |Wilson| ( |19^ . 


3.6. Solid [i2ci60]/[i3Ci®0] toward IRS 51 

The optical depth for the ice toward the fore¬ 

ground cloud traced by IRS 51 provides the only solid 
j12c16o]/[13c16o] ratio (68 ± 10) that can be com¬ 
pared to a gas-phase ratio along the same line-of-sight 
in this study. The ice-phase ratio is significantly lower 
than our derived gas-phase ratio of 165 ± 10. The 
[12 c160]/[13c16o] value in ice for IRS 51 is roughly 


consistent with other publ ished 


ues for ice in YSO targets: fBoogert et al. 

ratio of 71 ± 15 for t 


a 


/n^O] val- 
(2002a) derive 
le ice observed 
toward the massive YSO NGC 7538 IRS 9, using Keck- 
NIRSPEC, and a lower ratio of 45 ± 15 was observed 
toward t he CRBR 2422.8-3423 ed ge-on disk using VLT- 
ISAAC ( Pontoppidan et aT]|2005 ). 

4. DISCUSSION 

4.1. Galactic [^^C^®0]/[^^C^®0] heterogeneity 

We find significant heterogeneity in gas-phase 
ji2)gi6Qj/p3ci6Qj toward our targets (Figures 10 and 
[lT| ), with most values falling above local ISM values of 
65 — 69, as determined by sub-m illimeter observations 

The large spread 
either heterogeneity 


(Wilson 
in P^C 


Ibni \ 


0 ]/[ 


Milam et al. 2005). 

I'Sc^Qj 


suggests 


in the carbon isotopic composition in our local solar 
neighborhood (e.g. Rgc~ 8 kpc), and/or one or more 
new fractionation mechanisms for CO operating in the 
dense envelopes of young stars. The distribution of CO 
isotopologue ratios derived from our CRIRES data do 
not match the general Galactic trend, which suggests 
that the carbon reservoirs of young stellar envelopes 
do not necessarily isotopically represent t heir parent 


clouds, as is often assumed in GCE models (jPrantzos et 


al. 1996). 


[12C160]/[13C^60] abundance ratios in YSOs and 


molecular clouds (Fig ure 
Lambert et al.l ( 19941) an^ 


a F« 


In ultraviole t abso rption, 
federman et al. (2003) found 


high [^^C^'’Oj/p'^C^‘’0] values toward diffuse regions in 
the Ophiuchus cloud, with values ranging from ~ 117 
to 167. These results have been attributed to selective 
CO photodissociation by external ultravio let radiation 
(|Lanibert et al.ll994| |Federman et al.|2003|) . We discuss 
connections witn C(J photodissociation in Section |4.3[ 
Unusually high [^^C^®0]/[^^C^®0] values were also de- 
tected in rovibrational absorption using Subaru (A/AA « 
23, 000) toward the massive young stars LkHalOl and 
MonR2 IRS3, with values of ~ 137 and 158, respectively 


(Goto et al. 2003). Goto et al. (2003) suggested that 


are 


not the first to observe unusual 


photon-trapping in the rotational lines of the most abun¬ 
dant isotope as a possible explanation for the apparently 
elevated carbon isotope ratios, although the potential ef¬ 
fect was not quantified in their study. 

4.2. Ice-gas partitioning of ^^C^®0 and ^^C^®0 

4.2.1. Observed relation between the ice-gas fraetion and the 
CO isotopologue ratios 

A possible mechanism for the observed CO fraction¬ 
ation is mass-dependent exchange between gas and ice. 
Compared to single-dish, beam-averaged millimeter mea¬ 
surements of CO gas-phase isotopologue ratios in molecu¬ 
lar clouds, rovibrational absorption through protostellar 
envelopes is likely to probe denser gas. Dense, cold gas is 
known to lead to efficien t CO freeze-out onto dust grains 
(Jprgensen et al. 2004). If the adsorption/desorption 
rates of CO are mass-dependent, isotopologue fraction¬ 
ation may occur. Further, if the heavier isotopologue 
(^^CO) is more tightly bound to the dust grains, we 
would expect a increase in the [^^CO]/[^^CO] ratio in 
the gas, which is consistent with the detection of large 
columns of CO ice along the observed lines-of-sight. Here 
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TABLE 7 

CO ICE PROFILE OPTICAL DEPTHS, COLUMN DENSITIES, AND ICE FRACTIONS FOR THE CRIRES 

TARGETS'^’*’ 


Target 

Component 

'^mc 

(10^® cm 

(lO^® cm-2) 

CO ice fraction ^ 

DoAr24E 

single 

< 0.02 

< 0.01 

1.10 ±0.005 

< 0.01 

DoAr24E 

single 

< 0.02 

< 0.01 

1.09 ±0.003 

< 0.01 

VV CrA N'* 

low 

< 0.05 

< 0.03 

1.41 ± 0.02 

< 0.02 

VV CrA 

high 

< 0.05 

< 0.03 

7.62 ±0.07 

< 0.004 

VV CrA S'i 

single 

< 0.05 

< 0.03 

1.18 ±0.004 

< 0.03 

HL Tau'i 

low 

< 0.08 

< 0.05 

1.07 ± 0.02 

< 0.05 

HL Taud 

high 

< 0.08 

< 0.05 

4.47 ± 0.06 

< 0.01 

LLN 19® 

single 

0.03 ±0.002 

0.02 ±0.001 

2.02 ±0.03 

0.01 ± 0.001 

RNO 91 

low 

0.25 ± 0.01 

0.17 ± 0.004 

0.69 ±0.07 

0.19 ±0.005 

RNO 91 

high 

0.25 ± 0.01 

0.17 ± 0.004 

2.36 ±0.03 

0.07 ± 0.002 

Reipurth 50 

low 

0.51 ± 0.003 

0.34 ± 0.002 

1.98 ±0.16 

0.15 ±0.01 

Reipurth 50 

high 

0.51 ± 0.003 

0.34 ± 0.002 

10.8 ±0.77 

0.03 ±0.002 

IRS 63 

low 

1.69 ±0.01 

1.12 ±0.01 

1.90 ±0.03 

0.37 ± 0.004 

IRS 63 

high 

1.69 ±0.01 

1.12 ±0.01 

5.74 ± 0.15 

0.16 ±0.004 

IRS 51 

single 

4.54 ±0.05 

3.01 ± 0.03 

2.03 ±0.04 

0.60 ±0.01 


to form CO 9 ( 

Mennella et al. 12006 Oba et al. 12010 

lop- 

polo et al. 

2 UJ 

Lip. I'urther, we would expect the column 


Note. — 

^ Optical depths for all targets excepting LLN 19 are d erived from the middle co mponent as 
determined by the phenomenological decompositions in [Pontoppidan et al.| | f2003[ |. 
Uncertainties are Icr. 

" [ 12 Cl 60 ]lce/(PCl 60 ]l,e + [^^ClSOlGa.- 

Upper limits on the CO ice column from spectra in this paper. 

® Optic al depth is from t he middle component values from the phenomenological decomposi¬ 
tions in |Thi et ah] | |2010[ | . 

we discuss whether ice-gas partitioning is a possible ex¬ 
planation for the observed high CO isotopologue ratios 
in the gas. 

The relative fraction of CO sequestered in ice, av¬ 
eraged along the line-of-sight, is plotted against the 

S 1!^®0]/[^^C^®0] abundance ratios in the gas (Figure 
. While we might expect to see increased values of the 
abundance ratio for lines-of-sight with a high ice/gas 
fraction if ^^CO is preferentially bound in the ice man¬ 
tles, we instead find significant scatter, much larger than 
the individual errors on each data point, with no sim¬ 
ple apparent relation between these quantities. A corre¬ 
lation is particularly absent between cold-gas reservoirs 
and CO ice fraction, although these two quantities are 
likely cospatial and trace the same total CO reservoir. 

However, the line-of-sight with the largest ice fraction 
(IRS 51) also has the highest CO isotopologue ratio in 
the gas. Therefore, while other factors are likely equally 
important, we cannot rule out that ice-gas partitioning 
could play a significant role in the observed high gas- 
phase pCieoj/pCi^O]. 

With respect to other factors affecting the gas-phase 
ratios, we note that high [^^C^® 0 ]/[^^C^® 0 ] ratios are 
also ob served in the warm gas in some lines-of-sight (Fig- 
red diamonds), which is unexpected if ice-gas 


11 


ure 

partitioning is the only fractionation mechanism, since 
environments warm enough to fully sublimate the CO 
ice should reset the gas isotopologue fraction back to the 
canonical value of ^ 65. In order to reconcile this with 
ice-gas partitioning, there must be a mechanism to ef¬ 
fectively remove CO from the ice during the cold phase 
of envelope evolution. Chemical pathways removing CO 
from the ice include the successive hydrogenation of CO 


the gas in order for the ice to be a significant reservoir 
that could increase the gas-phase [^^C^®0]/[^^C^®0]. For 
IRS 51, the large discrepancy in [^^C^®0]/[^^C^®0] be¬ 
tween the ice (^ 68 ) and the gas (^ 165), along with a 
seemingly low CO ice fraction (~ 0.60), suggests instead 
that open-system chemistry may be occurring via CO ice 
reactions. 

While the data are inconclusive with regard to carbon 
isotopic fractionation between reservoirs of CO ice and 
gas, potential trends in cold- and warm-gas reservoirs can 
be further investigated with high-resolution spectroscopy 
for a larger sample of protostellar lines-of-sight. Further, 
models considering chemical evolution in protoplanetary 
environments including YSOs could consider interactions 
between ice and gas reservoirs as a potential influence on 
isotope partitioning. 

4.2.2. Comparison to theoretical and experimental studies 

To estimate whether there is an expectation of ice- 
gas partitioning for CO, we briefly review the theoreti¬ 
cal basis for a potential mass-dependence in the ice-gas 
balance. Follo wing a standard approac h for ice-gas par¬ 
titioning (e.g. |Tielens & Hagen |1982D, CO separations 
between the solid- and gas-phase reservoirs can be in¬ 
vestigated with respect to mass-dependent kinetics and 
the balance between the rates of ^^C ^®0 and ^^C ^®0 gas 
adsorption and desorption, as given by the relation. 


dn;. 


dt 


— .^ads Rde 


( 2 ) 


daka 


to formaldehyde (H 2 CO) and methanol (CH 3 OH) ( Hi- 
" ■ —- et 


:a et al.||2004[ IWatanabe fc Kouchi||2002[ ICuppen 
2009[ Fuchs et al. 2009p , or reactions with OH or 


where nice is the number density of ice, i?des = 
z/oexp(-Ag/fcr) x rtco,ice x (3, and Rads = nco,gas x 
''T-dust X Tid?' X i/S/cT/mco X / are the desorption and 
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Fig. 8 .— Rotational excitation diagrams for targets with little to no ice, corresponding to targets shown in Figures and[^. E 

is the energy of the Jth rotational state above the ground rotational state and k is the Boltzmann constant. One- and wo^emperature 
simultaneous fits were used as warranted by the rotational excitation structure. Derived high and low temperatures for each isotopologue 
are indicated on the plots. Error bars are la propagated from the Gaussian fits to the line profiles. Faded symbols are outliers that have 
been excluded from the fits. 
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adsorption rates, respectively. The factor, AH, is the 
binding energy of CO. The factor /3 accounts for CO 
desorbing only from the top monolayer of the ice, d is 
the grain radius, and / is the sticking coefficient. The 
mass of CO ( toco ) is an obvious mass-dependency of this 
relation, and vq is the pre-exponential factor. The num¬ 
ber densities of the CO gas and du st particles are given 
by nco.gas and Udust, respectively (Pontoppidan 2006). 
In addition to the mass itself, the binding energy and 
pre-exponential factors are likely both mass-dependent, 
although to unknown degrees. 

Some experimental work exists for the mass depen¬ 
dence of thermal desorption of CO. In an ultra-high vac¬ 
uum, te mperature-programmed desorption (TPD) ex¬ 
periment, Acharyya et al. (2007) measured binding ener¬ 
gies and pre-exponentials for desorption from pure ^^CO 
and ^^CO ices. The parameters were derived using an 
empirical kinetic model. However, no difference between 
the ^^CO and ^^CO binding energies could be measured 
within the experimental uncertainties: 858 ± 15 K and 
854 ± 10 K for and respectively. 

To estimate the degree of ice-gas fractionation allowed 
by the boundaries imposed by the experimental evidence, 
we solved for the ice abundance in Eq. assuming a 
10 K difference in binding energy between the two CO 
isotopologues and a pr e-exponential that scales as 
{Hasegawa et al. 1992). In Figuref^it is seen that, under 
these assumptions, the [^^CO]/[^^C] gas ratio may reach 
values as high as 170, but only for a narrow range of tem¬ 
peratures close to those of CO desorption (above which 
no ice exists). Essentially, a dust temperature exists at 
which ^^CO can efficiently remain bound to the dust sur¬ 
face, while ^^CO cannot. Nearly all of this effect is due 
to the assumed difference of only 10 K in the binding en¬ 
ergies. Compared to what is currently available, signif¬ 
icantly more precise laboratory measurements of mass- 
dependent binding energies would be needed to rule out 
binding-energy differences as a significant cause of ice-gas 
fractionation. While such temperature-dependent frac¬ 
tionation may be consistent with the observed spread in 
isotopologue ratios, it remains unclear to which degree 
the physical conditions of real molecular clouds and en¬ 
velopes must be fine-tuned to lead to significant ice-gas 
fractionation. 

More realistic modeling of ice-gas fractionation of CO 
is potentially complicated by several factors. For exam¬ 
ple, investigations of the effects of CO entrapment in 
the ice on the desorption and adsorption (both poten¬ 
tial mass-dependent processes) in astrophysical environ¬ 
ments have shown that diffusion of CO into the porous 
H 2 O ice leads to competing CO entrapment and diffusion 
effects on the desorption of C O, the degree to whi ch de¬ 
pends on the ice temperature ( Collings et al.|2003 1. Fur¬ 
ther, a recent study has shown that three-phase modeling 
of H 20 :C 02 ice desorption depends on the ice thickness 
and mixing rat ios, which may hav e significant effects on 
ice desorption (Fayolle et al.|20lT|. 


4.3. Oxygen isotopologues and CO self-shielding 

Isotope-selective photodissociation of CO due to self¬ 
shielding refers to the differential shielding of CO iso¬ 
topologues from photodissociation by far ultraviolet 
(FUV) radiation in proportion to their column densi¬ 


ties. With regard to oxygen isotopes, this leads to spa¬ 
tial regions characterized by more rapid destruction of 
^^C^^O and ^^C^®0 compared to the much more abun¬ 
dant ^^C^®0 molecule. Oxygen liberated during this 
process is thought to eventually end up in H 2 O, provid¬ 
ing an explanation for the overabundance of and 

H2^0 relative to in the early solar system. Isotope- 

selective photodissociation of CO due to self-shielding 
has been sho wn to be a relevant process in the interstel¬ 
lar m edium (van Dishoeck & Black 1988 Visser et al. 
200 ^. 


The solar system oxygen isotope anomaly is char¬ 
acterized by the mass-independent correlation between 
[^®0]/[^®0] and [^®0]/[^^0] found in meteoritic materials 
as compared to the mass-dependent distribution known 
in terrestrial rocks. Mass independence refers to rel¬ 
ative differences in [^®0]/[^®0] and [^®0]/[^^0] among 
primitive rocky objects that are nearly identical, sug¬ 
gesting changes in [^®0] relative to both [^^O] and [^®0], 
rather than the expected mass-dependent trend in which 
relative changes in [^®0]/[^^0] are about half those in 
[^®0]/[^*0] on a log scale. The solar system oxygen 
anomal y has defied conclusiv e explanation since its dis¬ 
covery ( Clayton et al. [1973 ), and has been investigated 
observationally by Brittain et al. (2005) and Smith et al. 


(20091, the latter study finding inass-mdependent signa- 


tures in oxygen isotopes observed in CO toward the VV 
CrA N disk, consistent with isotope-specific photodisso¬ 
ciation by CO self-shielding. 

The compiled oxygen data for this study, which also 
inclu des refined analyse s of VV CrA N and Reipurth 50 
from|Smith et al. (20091, are shown in Tableand Figure 
13 i'he mass-dependent lines in Figure [T^ are derived 
'Ey estimating the oxygen isotope ratios for the ISM as 
a function of distance from the Galactic center. These 
values can be estimated from the 4 kpc and 8 kpc values 
( Wilson|[l999 ), such that for [^®0]/[^®0]: 


[i"0]/[i®0]^^^ = (57.5 ± 10) Rgc + 97, (3) 

where the Icr uncertainty in the slope reflects the un- 
certainty in t he standard errors reported in Table 4 of 


Wilson (1999). For [^®0]/[^^0], equation § is combined 


wit h values for [^ ®0]/[^^ 0] in the ISM, which va ry from 
3.5 (|Wilson|1999|) to 4.1 (|Wouterloot et al.|2008|) (Figure 


We find a mass-independent signature in the oxygen 
isotopes relative to solar system abundances, consistent 
with CO self-shielding, toward the HL Tau disk, and the 
same mass-independent signature is retriev ed in the re¬ 
fined r esults for VV CrA N (as compared to Smith et al. 
(2009|)). These findings lend further support t or mod^ 
considering CO self-shielding on disk surfaces (iLyons & 


Young||2005t |Young et al.||2008|) . It is particularly inter¬ 
esting to note that the recent suggestions of planet for- 
mation in the HL Ta u disk based on new ALMA imaging 
( Tamayo et al.|20f5 ) would imply that CO self-shielding 
may be a process that occurs and/or persists later in the 
evolution of YSOs than previously thought. 

We would not expect UV photolysis to be a significant 
process in the foreground cloud traced by VV CrA S, or 
in the relatively embedded objects, RNO 91, Reipurth 
50, and IRS 63. Nevertheless, due to the errors in the 
rare isotopologue, a mass-dependent effect can- 
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Fig. 9. — Rotational excitation diagrams for targets with moderate to deep CO ice features, corresponding to targets in Figuresand 
E is the energy of the Jth rotational state above the ground rotational state and k is the Boltzmann constant. One- and two-temperature 
simultaneous fits were used as warranted by the rotational excitation structure. Derived high and low temperatures for each isotopologue 
are indicated on the plots. Error bars are la propagated from the Gaussian fits to the line profiles. Faded symbols are outliers that have 
been excluded from the fits. 
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Fig. 10.— Compilation of ratios vs. Rqc (kpc). Gray 

symbols show values from the literature, derived from dif¬ 

ferent carbon-bearing molecules, CO, CN, CH"^ and H 2 CO. Data 
from our study are labeled for YSO targets (filled circles) and 
foreground clouds (unfilled circles). Colo rs connect these data by 
molecular cloud with those in Figure [T^ For targets where there 
is both warm and cold gas along a line-ot-site, ratios are shown for 
the more robustly de termined warm gas, unless otherwise noted. 
Local ISM (68 di 15; iMilam et al. | (|2005|l) and solar system val¬ 
ues (86.8 di 3.8; [Scott et ai.| ( |2UU6| l) are noted. The black arrows 
indicate high ^^C/^'^C in diffuse r egions of Ophiuchus ||Lambert| 
jet al.||1994| Federman et aini2003 1 . YS O data for MonK,2 
LkH alUi/At C j.L 4ij U jcjoto et ai. 2003[ l and HL Tau (lower value; 
|Brittain et al.| l|2005|l) are marked (dark purple). 



study. For targets where both warm and cold gas was observed, 
CO ice fractions are plotted with red diamonds when warm gas 
was used in computing the fraction, and blue stars when cold gas 
was used. For targets where only cold gas was observed, CO ice 
fractions are shown with over-plotted blue ovals. Errors are Icr 
propagated from errors in ice and gas optical depths measurements. 
3(7 upper limits (down-facing arrows) are shown for DoAr24E N,S, 
VV CrA N,S, and HL Tau. 



Fig. 12.— The gas-phase [^^CO]/[^^CO] ratio predicted by the 
simple model described by Eq. assuming a difference in binding 
energy of 10 K between the two isotopologues, as allowed by the 
experimental uncertainty, and a pre-exponential that scales with 
The model is a snapshot at 5 X 10® yr and the individual 
curves are created by varying the gas density from 10® to 10® cm“® 
and temperature from 9.1 to 10.0 K, as indicated. 


not be distinguished from mass-independence in oxygen 
fractionation for RNO 91 and IRS 63. Additional high- 
quality CO data for YSOs would help in evaluating mass- 
independence in oxygen for objects in various stages of 
evolution. 

With regard to the rare CO isotopolgues, with the 
exception of low-temperature gas in Reipurth 50, the 
[®^C^®0]/[^^C^^O] ratios are within 1 — 2 (t of molec¬ 
ular cloud radio measurements for local ISM, further 
strengthening the claim that the solar system is anoma¬ 
lous i n ^^O/^'O with resp ect to the local solar neighbor¬ 
hood (Young et al. 2011). This consistency in our data 
is observational support for ®^C^®0 and ®^C^^O being 
similarly affected in young stellar environments. 


4.3.1. Evaluation of CO self-shielding on carbon isotope 

ratios 

While critical evaluation of mass-independence re¬ 
quires three isotopes, we qualitatively considered if 
isotope-selective photodissociation by CO self-shielding 
could be affecting the derived [®^C^®0]/[^^C^®0] toward 
our targets by comparing the oxygen and carbon isotope 
ratios where this was possible: VV CrA N, VV CrA S, 
HL Tau, RNO 91, Reipurth 50, and IRS 63. These are 
plotted in Figure [T^ with the qualitative relative trend 
for CO self-shielding overplotted on the data. The de¬ 
rived [^®0]/[^^0] changes by nearly a factor of two rel¬ 
ative to [i2Ci®0]/[i3ci60], while the [i^Ci^Oj/pC^sO] 
varies by roughly 1.5 times with [^^C]/[^^C]. If CO self¬ 
shielding were the dominant mechanism, we would ex¬ 
pect at least positive, direct changes in the carbon with 
the oxygen isotopes along an integrated line-of-sight, as 
well as similar variability in [®^C^®0] and [^^C^^O] com¬ 
pared with [i2ci6Qj/^i3ci 6Qj similar rates 

of photodissociatio n (e.g. van Dishoeck fc Black| 1988 


Visser et al. 2009|). Therefore, while we cannot be cer- 


tam that photodissociation of CO has no affect on the 
[®^C^®0]/[^^C^®0], we cannot claim that it is a dominant 
process affecting this carbon reservoir. 

Models for low-mass YSOs predict that 
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Fig. 13. — Comparison of oxygen isotope ratios between the in¬ 
terstellar medium and CO surrounding the targets in our study 
where 3 oxygen isotopes were observed. Oxygen isotope ratios 
are plotted in linearized delta notation, such that the [^^0]/[^^0] 
and [^®0]/[^®0] from local ISM values are shown in per mil us¬ 
ing the linearized form of the delta notation. Here, (5^^0' = 
103ln(i8Ri/i8_RLo,,i ism) and = [Cl80]/[Ci®0]. Ratios for 
the YSOs are derived from the best-constrained (high- or single-T) 
values from the rotational analyses (Tablej^. The solid and dashed 
lines show mass-independent fractionation (slope ~ 1) lines with 
different assumptions for the ratio, as indicated. The 

mass-dependent line (slope ~ 1/2) is also shown for comparison. 
Ellipses represent 65% confidence limits (standard deviation) de¬ 
rived from Icr errors from the Tq for each l ine. Targets are colored 
according to parent cloud (see also Figure 


should n ot increase significantly 


due to photodissociation (e.g. 

Visser et aL||2009 

Woods 

& Willacy 2009). Atomic 

carbon liberated 

by CO 


C’*' reacts rapidly to exchange carbon with CO by the 
reaction ^ + 35 K dWarin 

et_aI.| |T9Ml [van Dishoeck fc Black|[T988{ [Langer et al. 


19841, driving ^'^C back into CO rnd i diminishing the 


signature of selective photodissociation. However, very 
high [^^C^®0]/[^^C^®0] ratios derived for diffuse clouds 
have nevertheless been attributed to CO self -shielding. 
Values of ^ 117, 125 (Federman et al. 2003), and 164 
( Lambert et al.|19'M l have been observed toward diffuse 
Ophiuchus clouds that are thought to be subjected 
to high-UV fields, and where CO photodissociation is 
thought to dominate over the carbon exchange reaction 
that would otherwise erase a photochemical signature. 

While the upper range in p^C^®0]/[^^C^®0] in our 
sample is consistent with the Ophiuchus data (Figure 


10 ), we do not find the expected mass-independent sig¬ 


natures in oxygen isotopes consistent with CO self¬ 
shielding. Further, as compared to diffuse clouds, we 
would not expect the molecular gas deep inside embed¬ 
ded YSO envelopes to be subjected to such high UV 
fields given their r elatively high value s of visual extinc¬ 
tion {Ay ^ 10 (e.g. Vuong et al.||2003 )). 


5. CONCLUSIONS 

Heterogenous [^^C^®0]/[^^C^®0] gas-phase abundance 
ratios, with several unusually high values, have been 



Fig. 14. — Carbon isotope ratios against oxygen isotopic 
abundances in targets where all isotopologues of CO were ob¬ 
served. The [^^C^®0]/[^^C^^0] ratio is shown in purple, and 
j 12 (^ 16 qj^|- 12 q 18 qj derived [^®0]/[^^0] changes by 

nearly a factor of two relative to [^^C^®0]/[^^C^®0], while the 
[^®0]/[^^0] varies by roughly 1.5 times with [^^C^®0]/[^^C^®0]. 
The arrows indicate the relative trends we would expect if 
isotopologue-specific photodissociation by CO self-shielding were 
the predominant process affecting the oxygen and carbon isotope 
ratios. 

found using infrared absorption spectroscopy at very 
high resolution toward a sample of YSOs in nearby star¬ 
forming clouds. A possible partial explanation for this 
trend may be CO isotopologue partitioning between the 
ice and gas reservoirs. The degree of partitioning and 
the preservation of the relation in warm regions where 
the CO ice has sublimated would require the action of 
a process that can remove CO from the ice phase with¬ 
out returning it to the gas. Additional high-quality data 
could help substantiate these findings. 

Signatures consistent with isotopologue-selective pho¬ 
todissociation by CO self-shielding in oxygen were ob¬ 
served in the oxygen isotopes toward the disks, HL Tau 
and VV CrA. However, our results do not support CO 
self-shielding as a dominant mechanism for explaining 
the trends in the gas-phase [^^C^®0]/[^^C^®0] hetero¬ 
geneity. 

We conclude that chemical evolution models for young 
stellar environments should consider the interaction be¬ 
tween CO ice and gas reservoirs as a potential fractiona¬ 
tion mechanism affecting the carbon isotope budget and 
influencing carbon isotopic heterogeneity in the gas. In¬ 
teractions between these reservoirs could also contribute 
to isotopic heterogeneity between solar-type YSOs and 
their parent clouds. Finally, isotopologue-selective pho¬ 
todissociation by CO self-shielding should be considered 
an important mechanism affecting oxygen isotopes in 
protoplanetary disks, potentially persisting in the later 
stages of planet formation. 
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